Spin correlations and phase diagram of the perturbed Kitaev model 
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We present a general classification of the perturbations to the Kitaev model on the basis of their 
effect on it's spin correlation functions. We derive a necessary and sufficient condition for the spin 
correlators to exhibit a long ranged power-law behavior in the presence of such perturbations. We 
substantiate our result by a study of the phase diagram of the Kitaev model augmented by a loop 
term and perturbed by an Ising term, within a RVB mean-field theory. We estimate the stability of 
the spin-liquid phase against such perturbations and show that this model exhibits both confinement- 
deconfinement transitions from spin liquid to antiferromagnetic/spin-chain/ferromagnetic phases as 
well as topological quantum phase transitions between gapped and gapless spin liquid phases. 

PACS numbers: 75.10.Jm, 73.43.Nq 



The spin-1/2 Kitaev model on a honeycomb lattice 
(Fig 1), a rare example of an exactly solvable 2D quantum 
spin model, is known to have several interesting features 
[1-4]. The Hamiltonian for the model is given by 
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where j and I denote the column and row indices of the 
lattice and a^-y-^ are Pauli matrices. The model has spin 
liquid ground states and exhibits spin fractionalization 
at all energy scales with extremely short ranged (near- 
est neighbour) spin-spin correlations [3]. The spectrum 
is gapless for | Ji — J2I < J3 < Ji + J2 [1] and gapped 
otherwise. The gapped phase has Abelian anyonic ex- 
citations while the gapless phase supports non-Abelian 
anyons which could provide topologically protected sub- 
spaces for fault tolerant quantum computation [1, 7]. 

There have been several proposals for experimentally 
realizing this model in systems of ultracold atoms and 
molecules trapped in optical lattices, quantum circuits [5] 
and in layered iridates [6]. Almost inevitably such real- 
izations will have contaminating perturbations; therefore 
the effect of such perturbation on the spin-liquid phases 
of the Kitaev model deserves theoretical attention. Re- 
cently, the effect of one such perturbation, namely a Zee- 
man field given by the Hamiltonian = ~h'Y^^ erf, on 
the spin correlations of the model has been studied [8]. 
It was found that such a perturbation (arbitrarily weak) 
qualitatively alters the nature of the spin correlators by 
rendering them long-ranged (in power- law sense). How- 
ever, such an analysis has not been extended to other 
relevant perturbations. In particular, there is no general 
theoretical criteria to classify the perturbations according 
to their effect on the spin correlations. Furthermore, the 
stability of the Kitaev spin-liquid in the presence of per- 
turbations and the possible transitions of this perturbed 
Kitaev model from the deconfined spin-liquid state to the 
confined spin-ordered states has not been studied. 



In this letter, we study both the above-mentioned is- 
sues. First, we provide an explicit classification of differ- 
ent short-range perturbations, in the gapless spin liquid 
phase, based on the nature of the spin correlators they 
induce. We derive the necessary and sufficient condition 
for the spin correlators to exhibit a long ranged power- 
law behavior in the presence of such perturbations and 
show that the induction of a long-range spin-correlation 
by a magnetic field, studied in Ref. [8], constitutes a spe- 
cific example of this general condition. Following this we 
take up a concrete example of the second class of pertur- 
bation which does not induce long-range correlations and 
study the Kitaev model in the presence of a loop term and 
the Ising-like perturbative term . We analyze this model 
using a RVB mean-field approach and demonstrate that 
the model exhibits transitions from a deconfined spin liq- 
uid to confined Ising ordered antiferromagnetic (AFM), 
ferromagnetic (FM) or spin-chain (SC) phases. We also 
show that the model supports two distinct spin liquid 
phases with gapped and gapless deconfined spinon exci- 
tations with a topological quantum phase transition sep- 
arating them. We chart out the phase diagram of this 
model and estimate the stability of the deconfined phases 
as a function of the strengths of the loop and Ising terms. 
Both the issues studied here are highly experimentally 
relevant and our work constitutes a significant extension 
of our understanding of the perturbed Kitaev model. 

We begin with the study of the fate of spin correla- 
tions of the Kitaev model under most general class of 
short ranged perturbations. The total Hamiltonian is 
then given hy H = Hk + Hp, where Hp is a perturbing 
Hamiltonian made up of spin operators. Examples of Hp 
include Hz or the most general bilinear spin-spin inter- 
action term given by Hb = ^apcrfa^ including the 
experimentally most relevant anisotropic Heisenberg in- 
teraction Hh = Y.(tj) {^xx(7f(T^ + >^yy(j'^(7j + Kyafa^). 

In the analysis of the effects of these perturbations a 
crucial role is played by the Z2 flux operators defined as 
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FIG. 1: (Color online) Schematic representation of the Kitaev 
model on a honeycomb lattice showings the different links x, y 
and z and the two sublattices A and B. The sites labeled 1..6 
and their spin configuration is a schematic representation for 
a classical configuration with (Wp) = 1 where the green, blue 
and red spins point in the x, y and z directions respectively. 




FIG. 2: (color online) Left Panel: The fiux fiipped on applica- 
tion of of where i is the site at the center (circled): (a) a = z, 
(b) a = y {c) a = X (see Fig. 1 for definition of the x, y, z 
bonds). Right Panel: A schematic representation of y type 
covering VyR (red plaquettes) and VyB (blue plaquettes). A 
product of Wp over red or blue plaquettes yields E^. Anal- 
ogous coverings for x and z type coverings may be obtained 
by rotating the above figure by ±60° 



Wp = erf crf^lCTf CT|cr| (Fig. 1) on every plaquette. Since 
[Hk^ Wp] = these fluxes are frozen in the unperturbed 
model and this fact is crucial to the integrability of the 
unperturbed Kitaev Model. For the ground state of i?_R-, 
Wp = 1 Vp. However, Hp does not commute with the 
WpS and thus the Z2 fluxes acquires dynamics in the per- 
turbed model. The effect of this dynamics is key to deter- 
mining the nature of the spin correlators of the perturbed 
Kitaev model. To this end we define 3 conserved oper- 
ators, EE Wjiia'^) = riper^e where a = x,y,z 
and the product over WpS in the second expression is 
taken over a subset of plaquettes T^c hi the following 
way. Colour all the plaquettes with two colours, red and 
blue, as shown in Fig. 2 (right panel). Then TyniTyB) 
denotes product of all WpS on the red(blue) plaquettes 
which gives S^. Analogous coverings in the other two di- 
rections yields and . The S" operators correspond 
to global TT spin-rotations about the ath axis and with 
periodic boundary conditions we have S^S*' = (— 
where N is the total number of sites. 

With these definitions, wc now state and prove the nec- 
essary and sufficient condition that the perturbation Hp 
does not change at least one dynamic spin-spin correla- 
tion (i.e. ((((T"(i)(7Q (0)))) ) from short ranged to a long 
ranged one (in power-law sense). The condition is: 

[l]",i/p]_=0 Va, (2) 

where [..]_(+) denotes the commutator (anti- 
commutator). // this condition is violated, then at 
least one component of the dynamical spin-spin corre- 
lation becomes a power law. To prove this we take the 
case of ((cr^(t)crQ(O))) for concreteness (without loss of 
generality). We note that in terms of the Majorana 
Fermions we have af = —ibfa [1-3]. Thus the operation 
of erf creates a c Fermion on the i*'' site as well as two 
quanta of Z2 flux on adjacent plaquettes which share an 
Qf-type bond. This is schematically illustrated in Fig. 



2 (left panel). Thus the operators cr" and erg create 
such fluxes centered around r and respectively. Since 
the perturbation docs not conserve flux, one can always 
repeatedly apply Hp and move the pair of flux centered 
at r to the pair centered at and then annihilate them. 
This occurs at order ~ r of the perturbation theory and 
one has a string of perturbation operators joining and 
r. Such a process has a leading term A'' = e~'^^'^^^^'^^ 
which gives rise to exponential decay of spin correlations 
with correlation length ^ ~ l/ln(l/A). 

The above discussion indicates that in order to get 
long-ranged correlations, the perturbation term must 
neutralize the two pairs of fluxes centered at r and 
locally and independently without forming the interme- 
diate string whose length scales with r. This requires 
that the application of erf changes Y.^ for at least two 
values of /3, e.g. cr,f]+ = and [E^crf]- = 0. 

Thus the perturbing Hamiltonian, Hp, must be able to 
change the corresponding S^s in order to neutralize the 
flux locally. This is possible when the generic condition 

P^i/p]+=0 (3) 

is fulfilled for the corresponding two pair of /3s. Eqs. 2 (3) 
represents the necessary and sufhcient conditions that a 
perturbing Hamiltonian Hp docs not (docs) change the 
nature of the spin correlators. 

To illustrate these conditions via specific examples, we 
consider Hp ~ H^ (studied in Rcf. [8]). It is easy to 
sec that [cr;', = and [cr^,i;^]_ = 0. Thus the 

Z2 fluxes are neutralized locally leading to long-ranged 
spin correlators. For a Zeeman term in the z direction, 
only the zz correlators become power-law while xx or 
yy correlators remain exponential as before. The above 
results hold for Hp = Hi, when a ^ /3. Next, we consider 
Hp = Hi,. The individual terms in H}, are of the form 
afcj'j, where i and j are nearest neighbors. If (zj) is in 
the a"^ direction, they commute with all Wp. Otherwise, 
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each spin operator flips two distinct pairs of WpS on the 
plaqucttcs which share the bond i,i + Ca and j,j + Ca- 
However, it is easy to check that contains either both 
or none of the flipped WpS and hence commutes with 
all the three S". Hh therefore fails to destroy the short 
ranged nature of the spin-spin correlations. 

Having established the general setting, let us now 
look at the spin-spin correlation in detail. A 
typical n^^ order term in the perturbation ex- 
pansion for the spin-spin correlation is, r„ = 

• • • E(.„,„) / • • • / dr„(r(<(r)/i- ^ (ri) 

/i-^,,(T2).../iz„j„(T-„)a^(0)))/n!, where hf^ stand for 

the individual terms in Hp. Here we have taken 
all operators to evolve in Euclidean time, 0{t) = 
^Hkt q^-Hkt Consider now the case when r is suf- 
ficiently far from the origin. At any finite order, n 
(with n being even) T„ may be written as r„ = 

• • • E(.„,„) / rfn • • • / dT„(r(K(r)/i^^^ (n) 
^"2% (^2)-] [• • • /iw„(T„)(T^(0)] ))/«!, where we have di- 

vidcd up the series using [• • • ] such that the fluxes are 
neutralized around r and respectively by the opera- 
tors belonging to each group within the square brack- 
ets. Once the fluxes are neutralized, the 5f Fermions 
are no longer important except for an overall constant. 
The c Fermions now determine the details of the cor- 
relation. In the gapless phase of the Kitaev model the 
c fermions have a Fermi-surface and gapless excitations 
across it. Thus the c-Fermion propagator is an n"^ order 
free Fermion propagator, which in (2 + 1) dimension, is 
given by Oc{P,r,T) ~ (r^ -|- r^) '-"''/^'''^•'^ where p is the 
number of cr operators occurring in /i° . Thus we expect 
that the connected spin-spin correlation function goes as: 

{{ant)aSm)^X-{r'^tY''''/'^'\ (4) 

where A is the coupling constant. Eq. 4 reproduces the 
results of Ref. [8] where p — 1 and n ~ 2. 

We note that the above results are valid for infinite 2D 
systems and our conclusion may change for finite systems 
such as nanotubes, i.e., a cylinder of infinite length and 
(finite) perimeter L. In this case it may be possible to 
annihilate the Z2 flux locally by going around the finite 
direction of the cylinder even for a perturbation which 
does not induce power-law correlations in an infinite sys- 
tem. This leads to a crossover of the behavior in the cor- 
relation function. For L ^ r, it is easier to construct the 
string joining and r (at the lower order of perturbation 
theory) than going around the cylinder and the infinite 
geometry results hold. However, for r ^ L, it is easier to 
form a string round the axis of the cylinder annihilation 
of the Z2 flux leading to to power-law correlators. Such 
a term occurs at order ~ of the perturbation theory. 

Next, we study the phase diagram of the Kitaev model 
augmented by a loop term Hl = —i^J2p^p per- 



turbed by an Ising-like perturbation Hj = A^^-^^ ^i^j 
for Ji.2,3 = 1 and within a RVB mean-field theory. To 
this end, we use our earlier transformation af — —ibfci 
to map the spin model Hk + + Hj to a Fermionic 
model Hp. The resultant Hamiltonian becomes 

Hp = + H ^^l^j^^is:. 

jGA a—x.y links z link 

j ,k^plaquette zlink 
j a— all links 

where the subscript j, k S plaquette indicates that the 
sum is over sites which belong to the A sublattice of a 
given plaquette as schematically shown in Fig. 1. Note 
that for A = 0, the operators ibjbj, , commute with the 
Hamiltonian and are therefore a constants of motion. In 
this limit, H is exactly solvable. When A is turned on, 
these operators acquire dynamics and their fluctuations 
are ultimately expected to conflne the spinous through a 
confinement-deconfinement transition. 

To make further progress, we introduce RVB type 
mean-fields [3, 4] on the sites (corresponding to spin 
ordering) and on links (corresponding to the emergent 
gauge fields) of the hexagonal lattice: {ibjCj) = (cr|) ~ 
Ai(2), {-ibjCfJ = Pa, {"ibjb^.J = Ja, and {iCjC^'J = 70a. 
Note that keeping in mind the bipartite nature of the 
hexagonal lattice and to allow for possible AFM phases, 
we have introduced two mean-fields Ai and A2 corre- 
sponding to the two sublattices shown in Fig. 1 . Decom- 
posing the quartic term in Hp in using these mean fields 
we have the quadratic mean-field Hamiltonian, which, in 
momentum space, is given by 

k 

(/3 - 2k^,/j'o + e^*=^ + e'('=i+^^)) bfb? 

+ {icfbl~ibfcl){/3^il + \)+(3.,e^'^^ 

+P^e^'^^)-c,bfc^-C2bl^c?+h.c 

+K^l - (1 + A)73 7o. + (1 + 3A)AiA2, (6) 

where Jo (1 -f- A7^), aJo = (1-|-A)7^, Jq = A70X, 
/3Jo = (1 + A) 702, ci(2) = (l + 3A)Ai(2), and the momen- 
tum k = kiii + k2e2 with the unit vectors ei = i + y/VS 
and 62 = 2y/\/3. 

We now minimize H^f numerically and obtain the 
mean-field phase diagram of the model as a function of 
A and k. Note that the mean field solution is exact at 
A = 0. This phase diagram is shown in the left panel 
of Fig. 3. In accordance with our expectation, we find 
that at large positive (negative) A, the ground state of 
the model is an Ising AFM (FM) which corresponds to 
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FIG. 3: (Color online) Left Panel: The mean-field phase 
diagram for the model. The blue dashed lines represent 
confinement-deconfinement transitions. The triple points oc- 
cur at {\uK*i) = (1.0f,1.07) and (Aj.^a) = (3,5.64) for A > 
and (A3, Kg) = (-1.01,1.88) for A < 0. Right Panel; Plot of 
the loop order parameter (solid blue line), the spinon gap 
(red dashed line) and the FM and the AFM order parameters 
(black dashed lines) as a function of A for k = 7. 



confined phase of spinons while at small A, the model 
exhibits a deconfined gapless phase DCl. The transi- 
tion between DCl and AFM , at low k and positive A, 
occurs via a intermediate SC phase, which corresponds 
to antiferromagnetic alignment of spins along chains in 
X direction of the hexagonal lattice with ferromagnetic 
arrangement of such chains in the y direction. For nega- 
tive A, there is a direct transition to the FM phase from 
DCl (for small k). At high enough values of k, we find 
another gapped deconfined phase DC2. The transition 
between DCl and DC2 is and second order within mean 
field theory and is an example of a topological quantum 
phase transition. The confinement-deconfinement tran- 
sitions at high K always occur from DC2 to AFM/FM 
phases. These transitions are predicted to be first or- 
der within mean-field theory. The phase diagram ex- 
hibits two triple points at {XI, kI) = (1.01,1.07) and 
(A2,K2) = (3,5.64) for A > 0. These represent meeting 
points of AFM, SC and DCl and AFM, DC2 and DCl 
phases respectively. For A < 0, there is one triple point 
(A*, 4) = (-1.01,1.88) where the FM, DCl, and DC2 
phases meet. We also note that our mean-field analysis 
also gives an estimate for the stability of the deconfined 
phase of the Kitaev model (-0.07 < Ac < 0.08 for k = 0) 
under external perturbing Ising term which may be im- 
portant for physical realization of the Kitaev model and 
for quantum computing proposals based on it [5, 7]. 

The plot of the loop order parameter (Wp), the spinon 
gap, and the AFM and the FM order parameters as ob- 
tained from the mean-field theory, is shown, for /t = 7, as 
a function of A in the right panel Fig. 3. We note that all 
the order parameters show discontinuous changes at the 
transition points indicating first order transitions. The 
spinon gap, in contrast, increases linearly and contin- 
ually with A indicating a second order quantum phase 
transition between DCl and DC2 phases. The presence 
of this topological quantum phase transition and the lin- 
ear variation of the spinon gap with A can be understood 
qualitatively from iJ,nf- For large k, it requires a large 



A to destabilize the Kitaev ground state in favor of Ising 
AFM/FM. In addition, numerically we find that in the 
Kitaev phase 72 (7^) ~ 1(0). As a result, beyond a crit- 
ical value of A = Ac, the effective couplings along the 
links, Ji_2 ~ (1 + Xjx), J3 ^ 7z(l + A), fail to satisfy 
\Ji — Jzl < J3 < Ji + J2 thus leading to a gapped phase 
via a topological quantum phase transition [1, 2]. The 
spinon gap in this gapped phase varies linearly with J3 
[2, 3] and hence shows a linear variation on A. At small k, 
the confinement-deconfinement transitions to the SC/FM 
phases occur before Ac is reached and hence the topolog- 
ical phase transitions do not occur. 

To conclude, we have presented a general classification 
of the contaminating interactions of the Kitaev model 
based on their effects on the spin correlators of the model 
and have derived a necessary and sufficient condition 
for the interaction to induce power-law spin-spin corre- 
lations. We have also presented the phase diagram of 
the Kitaev model, augmented by a loop term and per- 
turbed by an Ising Hamiltonian, and have shown that the 
model exhibits a rich phase diagram with several inter- 
esting transitions. Our estimate suggests that the topo- 
logical phase of the Kitaev model is unstable to about 
10% contamination by Ising interactions. 
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